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Abstract 

The production of wrong sign charmed mesons b — > D^X, £)(„) = (D°, D + , D s ), 
is studied using the data collected by the DELPHI experiment in the years 1994 
and 1995. 

Charmed mesons in Z — > bb events are exclusively reconstructed by searching 
for the decays D° -> K~ir + , D + -»■ K-n + n + and D+ -> 0tt+ -> K + R-ix + . 
The wrong sign contribution is extracted by using two discriminant variables: 
the charge of the 6-quark at decay time, estimated from the charges of identified 
particles, and the momentum of the charmed meson in the rest frame of the 
6-hadron. 

The inclusive branching fractions of 6-hadrons into wrong sign charm mesons 
are measured to be: 

B(b -> D°X) + B{b -> D~X) = (9.3 ± 1.7 (stat) ± 1.3(syst) ± 0A(B))% , 
B(b^D-X) = (10.1±1.0(stat)±0.6{syst)±2.8{B))% 

where the first error is statistical, the second and third errors are systematic. 
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1 Introduction 

Decays b — > c are expected to occur through the Cabibbo favoured transitions b — > 
cVr _ and W~ — > cs 1 . Hence, neglecting b — > -u transition and .D mixing, 6-hadron 
decays to wrong sign charmed mesons are in fact double charm transitions. The double 
charm rate is related to n c , the mean number of charm quarks (and anti-quarks) produced 
per 6-decay: 

n c = 1 — Bib — > no open charm) + 2B(b — > charmonium) + B(b — > double charm) (1) 

which can be predicted by Heavy Quark Effective Theory (HQET) based calculations of 
the semileptonic B meson branching fraction [1]. 

Evidence for wrong sign charm production and double charmed 6-decays has been 
found both at the T(4S) and at LEP. ARGUS [2] and CLEO [3] have shown evidence 
for the two-body transitions B — > D^ + D^ 2 . From the analysis of the D s momentum 
spectrum, these decays are found to contribute about half of the total D s production at 
the T(4S f ), the remainder coming from either B -> D ( * )+ D** or B -> D { s * )+ D^ 7T, p, UJ 
(where D** denotes an orbitally excited D meson). By using D-lepton correlations, 
CLEO has observed wrong sign D production [4] . ALEPH has reported evidence for b — > 
DD( S )X decays with both charmed mesons reconstructed [5]. The observed DDX signal 
is shown to originate either from B — > D (*)_dW K W or from B — > D**D with a subsequent 
decay of the orbitally excited state D** into D^K. 

In this paper, the DELPHI data are used to measure the inclusive branching fractions 
of 6-hadrons into wrong sign charm mesons, B(b — > DX) and B(b — > D~X). D°, D + and 
Df mesons are exclusively reconstructed in Z — > 66 events, recorded by DELPHI in the 
years 1994 and 1995. The wrong sign contribution is extracted by using two discriminant 
variables: the charge of the 6-quark at decay time, estimated from the charges of identified 
particles, and the momentum of the charmed meson in the rest frame of the 6-hadron. 

2 Experimental procedure 

2.1 The DELPHI detector 

A detailed description of the DELPHI detector and its performance can be found in 
reference [6]. Only the subdetectors relevant to the present analysis are described in the 
following. 

The tracking of charged particles in the barrel region is accomplished with a set of 
cylindrical tracking detectors whose axes are oriented along the 1.23 T magnetic field and 
the direction of the beam. 

The Time Projection Chamber (TPC), the main tracking device, is a cylinder of 30 cm 
inner radius, 122 cm outer radius and a length of 2.7 m. For polar angles between 39° 
and 141°, it provides up to 16 space points along the charged particle trajectory 3 . 

The Vertex Detector (VD), located nearest to the LEP interaction region, consists of 
three concentric layers of silicon microstrip detectors at average radii of 6.3 cm, 9.0 cm 
and 10.9 cm. Since 1994, the innermost and the outermost layers were equipped with 
double sided silicon microstrip modules allowing both R<p and z measurements. 

1 Charge conjugate reactions arc implied throughout this paper. 

2 In the following, D (0( s )) denotes cither D° or D+ (D°, D+ or Df). 

3 In the DELPHI frame, the z axis is defined along the electron beam direction, the x axis points towards the centre 
of the LEP ring and the y axis points upwards. The polar angle to the z axis is called 6; the azimuthal angle around the 
z axis is referred to as <j>. The radial coordinate is R = \J x 1 + y 2 . 



2 



Hadrons are identified using the specific ionization (dE/dx) measured in the TPC and 
the Cherenkov radiation detected in the barrel Ring Imaging CHerenkov counter (RICH) 
placed between the TPC and the Outer Detector (OD). 

2.2 Event sample 

For this analysis, the data collected by the DELPHI experiment in the years 1994 and 
1995 at ^/s close to 91.2 GeV are used, corresponding to about 2.1 million hadronic Z de- 
cays. Simulated hadronic events are generated with the JETSET 7.3 program [7]. Full 
detector simulation is applied to Monte Carlo events which are subsequently processed 
through the same analysis chain as the real data [6]. 

The decays B -> D^+D^, B -> D { s * )+ D**, B -> D { * )+ D^it, p, w, B -> D**D -> 
D&KD and B -> £)(*)£)(*) are used to model 6-decay into wrong sign charmed 
mesons. The 6-hadron decay properties to right sign charm are adjusted to match the 
latest experimental values [8]. In total, a sample of about 58,000 b — > D^X and about 
99,000 b — > D( S )X events, with forced to decay into the modes used in the analysis, 
has been generated. The background is modelled with about 3.2 million Z — > qq and 
about 1.8 million Z — > bb Monte Carlo events. 

Hadronic Z decays are selected by requiring at least five charged particles and a total 
charged energy larger than 12% of the collision energy [6]. The tagging of 6-quark jets is 
based on four discriminant variables, the most important one being the probability for all 
tracks to originate from the primary interaction vertex, calculated from the track impact 
parameters with respect to this point [9]. The other variables are defined for jets with 
a secondary vertex: effective mass of the system of particles attached to the secondary 
vertex, rapidity of these particles with respect to the jet direction and fraction of the 
charged energy of the jet included in the secondary vertex. All jet ft-tags in the event 
are combined and the cut on the event probability is chosen such that about 90% of 
the reconstructed charmed mesons originate from 6-hadron decay. Correspondingly, the 
Z — > bb selection efficiency varies between 58% and 74% for the different charm modes. 

Each selected event is divided into two hemispheres by the plane orthogonal to the 
axis of the most energetic jet and passing through the primary interaction point. 

2.3 Charmed meson reconstruction 

Charged particles are selected as follows: momentum larger than 100 MeV/c, relative 
error on the momentum measurement smaller than 100% and Rip (z) impact parameter 
with respect to the primary interaction vertex smaller than 4 cm (4 cm/sin 

Charmed mesons are searched for in the decay modes D° — > K~tt + , D + — > K~h + 'k + 
and Df — > <fiir + — > K + K~n + by trying all possible combinations of charged particles in 
the hemisphere. The dE/dx values of the kaon and pion candidates are required to be 
consistent with the respective mass hypotheses. For D + K~7t + tt + decays which suffer 
from a high level of combinatorial background, the kaon must be tagged additionally by 
the RICH. To allow for a precise reconstruction of the decay vertex, at least two 
tracks in each combination are required to have associated hits in the Vertex Detector. 

Track combinations satisfying these criteria are fitted to a common vertex. The \ 2 ' 
probability of the fit must exceed 0.01%. Combinations containing a fragmentation track 
are rejected by requiring the D (D s ) vertex to lie at least three (two) standard deviations 
away from the primary interaction point and imposing the requirement Xe > 0.15 on the 
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energy fraction x E = E D(g) / E beam . For Df — > 07r + , a selection at ±12 MeV/c around 
the nominal mass is applied to the reconstructed K + K~ mass. 

D candidates are selected by using four discriminant variables: the RICH information 
for the kaon candidate, the decay length from the primary to the charm vertex divided by 
its error, the energy fraction xe and the cosine of the charm decay angle 9d, defined as the 
angle between the K momentum vector in the D meson rest frame and the D momentum 
vector in the laboratory frame. The cos#d distribution is flat for the signal and peaked 
at —1 for the combinatorial background. For D s candidates, two additional variables are 
used: the reconstructed K + K~ mass and the cosine of the helicity angle 9h- The latter 
is defined as the angle between the K + and the D s direction in the rest frame. The 
signal follows a cos 2 9h distribution while the background is flat in cos 9h- The different 
variables %i £1X6 combined by using a likelihood ratio: 

Yin \ R ( D (*)) mn \ TT S ^ D (°)) (0 \ 
X{Dw) ~ T+WD£) > R ^-ii Bl ( Xi ;D {s) ) (2) 

where Si and Bi are the normalised distributions of Xi for the signal and the combinatorial 
background, respectively, as obtained from the simulation. The combined variable is 
defined such that the target value is X = 1 for the signal and X = for the background. 
For each decay mode, the selection cut on the variable X is adjusted on simulated events 
to optimise the statistical significance of the signal. The following selections are found: 
X(D°) > 0.8, X(D+) > 0.6 and X(D S ) > 0.95. 

For each selected candidate, the invariant mass of the track combination is computed 
(Figure 1). 7345 (6906, 984) D° (D + , D s ) candidates are found in the signal window 
corresponding to an interval of about ±2<r around the signal peak. The remaining com- 
binatorial background is determined by a fit to the invariant mass distribution. The 
fit uses a Gaussian function for the signal and a linear parametrisation for the com- 
binatorial background. The satellite peak, due to the D + — > K + K~n + decay, in the 
D+ — > 07r + channel is also fitted by a Gaussian function. In this way, the combina- 
torial background is found to be 3038 ± 43, 4677 ± 66 and 404 ± 12 for D° -> K~Tr + , 
D + — > K~n + 'K + and — > 07r + , respectively. 



2.4 Discriminant variables 

The discriminant variables used for selecting wrong sign decays are constructed by 
using a common DELPHI analysis package called BSAURUS. Details on how the different 
BSAURUS variables are formed can be found in reference [10]. 

The flavour of the D( s ) meson, i.e., the charge of the c-quark confined in the charmed 
meson, is determined from the charge of the kaon for the channels D° — > K~n + and 
D + — > K~7T + 7T + , and from the charge of the pion for the channel — > 07T + . 

The charge of the 6-quark at decay time in the hemisphere of the charmed meson is ob- 
tained from the BSAURUS Decay Flavour Neural Network (BDFNN). The approach is to 
first form the weighted sum of particle charges in the hemisphere excluding the particles 
from the exclusive decay of the , in order to avoid a possible bias. The weighting factor 
is constructed from the conditional probability for the track to have the same charge as 
the decaying 6-quark, and is determined via a neural network technique based mainly on 
particle identification variables for kaons, protons, electrons and muons and B — D vertex 
separation variables. In order to use optimally the event information, the resulting hemi- 
sphere charges are constructed separately to estimate the 6-charge at both production 
and decay time and this is repeated for each of the 6-hadron types (B + , B°,B®, 6-baryon). 
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In a final step, these hemisphere charges form the main input variables to a neural net- 
work trained to find the 6-quark charge in combination with BSAURUS 5-hadron type 
tagging probabilities and also including charge correlation information from the opposite 
hemisphere. 

The wrong sign tag Y(D( S )), the first discriminant variable, is obtained by correlating 
the BDFNN output with the flavour of the charmed meson: 

f BDFNN iorD (s) 
Y W*)) ~ \ 1- BDFNN for D {s) ' [6) 

The target values of the BSAURUS Decay Flavour Neural Network are BDFNN = 1 and 
BDFNN = for b and 6-hadrons, respectively. Hence, the target value of the wrong sign 
tag is Y = for wrong sign decays and Y — 1 for right sign. 

While wrong sign and B s right sign decays contribute about equally to the D s pro- 
duction in Z — > bb events, the wrong sign production mechanism is strongly suppressed 
in the case of D mesons. Hence, for selecting wrong sign D° and D + mesons, stronger 
discrimination is required and a second variable, the momentum of the D meson in the 
6-hadron rest frame, p(D), is used. 

The 6-hadron four-momentum in the hemisphere of the charmed meson is inclusively 
reconstructed in BSAURUS using the following procedure. An initial estimate of the 
6-hadron momentum p raw , energy E raw and mass m raw = a/ E? aw — vlaw ^ s formed from 
particles with high rapidity for events with more than two-jets and from the sum of "&- 
weighted" four-vectors for the two-jet case. This weighting involves the use of neural 
networks trained to identify tracks originating from the weakly decaying 6-hadron in 
the hemisphere. E raw is then corrected, hemisphere-by-hemisphere, motivated by the 
observation (in Monte Carlo simulation) of a correlation between the energy residuals 
AE = E raw — E true and m raw (which is approximately linear in m raw ) and a further 
correlation between AE and Xh = E hem / E beam , where E hem is the sum of the energies 
of all particles reconstructed in the hemisphere, resulting from neutral energy losses and 
inefficiencies. These effects are parametrised and corrected for, after which the resolution 
obtained in p(D) is about ±300 MeV/c. 

The two discriminant variables are shown in Figure 2. 

2.5 The fit 

For each charm decay mode, the numbers of wrong sign and right sign events, Nyy 
and Nr, are determined by a fit to the above-mentioned discriminant variables. The fol- 
lowing components can contribute to the distributions of these variables in the real data: 
wrong sign b — > D^X mesons, right sign b — > D^X mesons, meson background 
(contamination by charmed mesons produced in Z — > cc events) and combinatorial back- 
ground. The shapes of the distributions of these four components (F w , F R , F cc and 
pBkgrd^ are determined from the Monte Carlo simulation. In each fit, the number of 
charmed mesons from Z — > cc events is fixed to the value calculated from the partial 
width R c = 0.1702 ± 0.0034, the fragmentation probabilities f(c -»• D°) = 0.552 ± 0.037, 
/(c -> D+) = 0.237 ± 0.016, f(c -> D s ) = 0.121 ± 0.025 [11] and the acceptance deter- 
mined from the simulation. The numbers are found to be 436 ± 30, 266 ± 19 and 73 ± 15 
for D° — > K~7i + , D + — > K~7t + tt + and — > 07r + , respectively. The normalisation of 
the combinatorial background is fixed to the values quoted in Section 2.3. 
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Sample 


Wrong sign evts. 


Right sign evts. 


£w/ e R 


B(b^D (s) X) (W s 
B(b^D (a) ,D {3) X) y /0 > 


D° -> K-it + 

D + -> ^-7T + 7T + 
£>+ -> 07T+ 


383 ± 81 
186 ± 86 
286 ± 42 


3, 396 ± 110 
1,811 ± 101 
221 ± 39 


0.92 ±0.02 
0.80 ±0.03 
1.01 ±0.03 


11.0 ±2.1 ± 1.5 
11.4 ±4.7 ±3.5 
56.2 ±5.7 ±3.3 



Table 1: The fitted numbers of wrong sign and right sign mesons, the relative selection 
efficiency of wrong sign and right sign mesons and the fraction of wrong sign events in the 
charm signal. The error on the number of events is purely statistical. The error quoted 
on ew/tR is just that due to Monte Carlo statistics. The first error on the wrong sign 
fraction is statistical; the second one is the sum of all systematic uncertainties listed in 
Table 2. 



Selected D s meson candidates are arranged in 10 bins, i, of Y(D a ) (bin width 0.1) and 
the resulting one-dimensional histogram is fitted by the function: 

N t = N W F™ + N R F t R + N C - C F? + N Bkgrd F* k9rd . (4) 

The normalisations SjF^ = 1, Sji^ = 1, SjFj CC = 1 and Y] i Ff k9rd = 1 are used. Selected 
D° and D + candidates are arranged in 4 bins, i, of Y(D) (bin width 0.25) and 13 bins, 
j, oip(D) (bin width 200 MeV/c) and the fit function: 

N tj = N W F% + N R F* + N C5 F% + N Bk9rd Ff grd (5) 

is used. The fit algorithm accounts for finite Monte Carlo statistics [12] and the total 
number of selected candidates is used as a constraint. By applying the algorithm to 
simulated Z — > qq events, no significant bias in the fit result is observed. 

The results obtained by fitting the real data are shown in Figures 3, 5 and 6. The 
Monte Carlo model of the combinatorial background is tested on real data D< s \ candidates 
selected outside the signal mass window (Figure 4). The numbers of wrong sign and right 
sign events for each decay channel are given in Table 1. For the one-dimensional fit, the 
value of the x 2 is 4.8 compared to 10 — 1 degrees of freedom. The two-dimensional fit has 
52 — 1 degrees of freedom and the x 2 is 52.0 (62.3) for the D° (D + ) sample. 



3 Results and systematic uncertainties 

For each charm decay mode, the fraction of wrong sign events b — > D^X in the 
signal b — > D^,D^X is calculated from the numbers of wrong sign and right sign 
events: 

B(b -> D {S) X) N w 
B(b -> D (s) , D {S) X) N w + {e w /e R )N R ' 1 j 

The results are given in Table 1. The factor ew/^R in the denominator of Eq. 6 corrects for 
the different selection efficiencies of wrong sign and right sign mesons and was obtained 
from the simulation. 

The wrong sign charm model used for the fit assumes a 50% contribution of two-body 
decays B — > D (*)+£(*) to the wrong sign D s signal and the same relative contribution 
of B — > D**D to wrong sign D. The corresponding modelling uncertainty (Table 2) is 
estimated by varying these ratios within (50 ± 13)% and (50 ± 25)%, respectively. 



Source 



Value 



^■B{b^D ,D°X) y /0 ' 



A 



B(b^D-X) 
B(b->D±X) 



(%) 



A 



Bjb^DJX) 
B(b->DfX) 



(%) 



Model dependence (w.s.) 

B -> D**D fraction 

B -> D ( s * )+ D^ fraction 
Model dependence (r.s.) 

B(6 -> D°l~vX) 

B(b -> D+l-VX) 

B(b -> D+1-9X) 

B(b -> D°D~X) 

B(b -> D+D-X) 

B{b -> D+DJX) 

B(b -> 

B(6 -> D + DX) 

B(b -> £>+M) 
Z — > cc background 
Combinatorial background 

normalisation 

shape 

Total 



(50 ±25)% 
(50 ± 13)% 



(6.60 ± 
(2.02 ± 
(0.87 ± 
(9.10 ± 
(4.00 ± 
(1.17 ± 
(6.45 ± 
(1.80 ± 
(1.17 ± 



0.60)% 
0.29)% 
0.28)% 
3.35)% 
2.05)% 
0.71)% 
2.08)% 
0.96)% 
0.71)% 



0.12 

0.20 

0.08 

1.40 

0.09 

0.37 

0.41 
1.53 



1.15 

0.61 

0.59 

2.18 

0.05 

1.46 
1.56 
1.04 

3.53 



0.54 

0.61 

2.09 

2.39 
0.56 

0.29 

0.13 
3.34 



Table 2: Breakdown of the systematic error on the wrong sign fractions. For the total, the different components have been added 
quadrature. 
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The different decays used to model the right sign component are collected into four 
categories: b — > D( s )l~&i(X), b — > D^ir,p, u, . . . , b — > D^Dj(X) and 6 — > D^D(X). 
To estimate the systematics related to the right sign modelling, the relative contributions 
of 6 — > D(_ 8 )l~i>i(X), 6 — > D( S )D~(X) and 6 — > D^D(X) to the signal are varied. The 
ranges are obtained from recent measurements (Table 2). The weight of each category is 
varied separately and, for the total, the different contributions are added in quadrature. 

Further contributions to the systematic error are: normalisation of the Z — > cc back- 
ground (uncertainty in R c and in the fragmentation probabilities), normalisation of the 
combinatorial background (uncertainty of the fit to the invariant mass distribution) and 
uncertainty in ew/tR- For the D + — > K~n + n + sample which is particularly affected 
by the combinatorial background, instead of using simulated data in the signal window, 
the fit is repeated using real data candidates selected outside the signal window. Both 
approaches are statistically consistent and the difference in the fit result is interpreted as 
a systematic uncertainty related to the combinatorial background shape. 

4 Conclusion 

The production of wrong sign charm mesons in 6-hadron decay, b — > D^X, = 
(D°, D + , D s ), was studied using the DELPHI data collected in 1994 and 1995, leading 
to a measurement of the fraction Bib — > Dt s \X)/B(b — > Di s \,Di s \X). Combining this 
measurement with the branching fractions B(b — > D°,D°X) = (60.5 ± 3.2)%, B(b — > 
D ± X) = (23.7±2.3)% and B(b -> DfX) = (18 ±5)% [8], the following result is obtained 
for wrong sign D: 

B{b -> D°X) + B{b -> D~X) = (9.3 ± 1.7(stat) ± 1.3(syst) ± 0A(B))% . 

The first uncertainty is statistical, the second one is the sum of all systematic errors listed 
in Table 2 (accounting for correlated model systematics) and the last one corresponds to 
the uncertainties in B(b — > D°, D°X) and B(b — > D ± X) (note that the quoted statistical 
error includes both real data and Monte Carlo statistics). This value is in good agreement 
with previous measurements by CLEO [4] and ALEPH [5]. The inclusive branching 
fraction for wrong sign D s is found to be: 

B(b -> D~X) = (10.1 ± l.O(stat) ± 0.6(syst) ± 2.8(B))% . 

Again, the first uncertainty is statistical, the second one is the total systematic error of 
Table 2 and the last one corresponds to the uncertainty on B(b — > DfX). This value 
agrees with the total D s production rate at the T(4S), B(B -> DfX) = (10.0±2.5)% [8], 
where the dominant source of D s production is double charm 6-decay. 
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Figure 1: The invariant mass of selected D° — > i^ _ 7r + , D + — > fT~7r + 7r + and D+ — > 07r + — > K + K~ir + candidates. The points with error 
bars are the real data. The solid line is the result of the fit mentioned in the text. The signal window is shown by dashed vertical lines 
and the dashed histogram represents the Monte Carlo expectation for the combinatorial background (arbitrary normalisation). 
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Figure 2: The wrong sign tag Y(D^) and the momentum of the charmed meson in the 
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Figure 4: Same as Figure 3 for candidates selected outside the signal mass window. 
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Figure 5: The D° momentum in the 6-hadron rest frame p(D°) in bins of the wrong sign 
tag Y(D°). The data are the points with error bars; the histograms are the components 
of the fit function (as described in the text). 
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Figure 6: Same as Figure 5 for D + . 



